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Assuming that cosmological dark matter consists of weakly interacting massive particles, we use the 
recent precise measurement of cosmological parameters to predict the guaranteed rates of produc- 
tion of such particles in association with photons at electron-positron colliders. Our approach is 
based on general physical principles such as detailed balancing and soft /coUinear factorization. It 
leads to predictions that are valid across a broad range of models containing WIMPs, including su- 
persymmetry, universal extra dimensions, and many others. We also discuss the discovery prospects 
for the predicted experimental signatures. 



Introduction — A variety of astrophysical and cosmo- 
logical observations indicate that a substantial fraction 
(perhaps as much as 30%) of the energy density in the 
Universe is due to non-relativistic, non-baryonic, non- 
luminous matter. The microscopic nature of this "dark" 
matter is at present unknown. Perhaps the most attrac- 
tive explanation is provided by the "WIMP hypothesis": 
dark matter is assumed to consist of hypothetical sta- 
ble particles with masses around the scale of electroweak 
symmetry breaking, in the 10 GeV - 1 TeV range, whose 
interactions with other elementary particles are of the 
strength and range similar to the familiar weak interac- 
tions of the Standard Model. Such weakly interacting 
massive particles (WIMPs) naturally have a relic abun- 
dance of the correct order of magnitude to account for 
the observed dark matter, making them appealing from 
a theoretical point of view. Moreover, many extensions 
of the Standard Model contain particles which can be 
identified as WIMP dark matter candidates. Examples 
include supersymmetric models Ij , models with "univer- 
sal" extra dimensions 0, little Higgs theories 0, etc. 

Particle accelerators of the next generation, such as the 
Large Hadron Collider (LHC) at CERN and the proposed 
linear electron-positron collider, may have enough energy 
to produce WIMPs. Once produced, WIMPs escape the 
detector without interactions, leading to an apparent en- 
ergy imbalance, or "missing energy" signature. In this 
letter, we will use the known cosmological abundance of 
WIMPs to predict the rate of such events at future col- 
liders. In sharp contrast to all existing studies, we will do 
so without making any assumptions about the details of 
the particle physics model responsible for the WIMP: our 
results are equally valid in all theories listed in the previ- 
ous paragraph, as well as in any other model containing 
WIMPs. In our approach, all model uncertainties are 
reduced to a single continuous parameter with a trans- 
parent physical meaning, plus a small number of discrete 
choices such as the spin of the WIMP. 

WIMP Annihilation Cross Sections from Cosmology — 
We assume that the observed dark matter is entirely due 
to a single WIMP particle, x- This particle carries a 
new conserved quantum number which prevents it from 



decaying into lighter Standard Model particles, making 
it stable. At the same time, two WIMPs can annihilate 
into a pair of Standard Model particles Q: 

x + x^x, + x,, (1) 

where Xi — l,q,g, . . . can be any Standard Model parti- 
cle. We assume that Q is the only process important for 
the determination of the x relic abundance; that is, we 
ignore the possibility of coannihilations between x's and 
other exotic particles. The present dark matter density 
depends on the cross sections of reactions in the limit 
when the colliding x particles are non-relativistic. If v is 
the relative velocity of two x's, each cross section can be 
expanded as 

oo 

a.. = 5:ai^).2^ (2) 

j=o 

where cr'^*'^ only receives a contribution from s-wave an- 
nihilation, cr*^^^ receives contributions from s- and p- wave 
channels, etc. It is clear that at low v, the lowest non- 
vanishing term in Eq. |(5Jl will dominate. We define 

aa„ = ^a(^»), (3) 

i 

where Jq is the angular momentum of the dominant par- 
tial wave contributing to x annihilation in a given model, 
and the sum only runs over the final states that give 
a non- vanishing contribution at this order in J. For 
most of this letter, we will restrict our analysis to two 
cases: Jq — Q and Jq = 1. We will refer to WIMPs 
in each case as "s-annihilators" and "p-annihilators" , re- 
spectively. (The extension of our analysis to higher Jq 
is trivial but rather unmotivated since all known models 
predict either s-wave or p-wave annihilation.) 

The present cosmological abundance of WIMPs is 
mainly determined by the values of Jq and (Tan, with 
only a weak dependence on other parameters such as the 
WIMP mass M^, its spin S*^, etc. It is therefore these 
parameters that are strongly constrained by cosmologi- 
cal observations. In Figure ^ we show the cosmological 
constraint on a^n, as a function of M^, for s- and p- 
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FIG. 1: Values of the quantity (Tan defined in ^ allowed at 
2a level as a function of the WIMP mass. The lower (up- 
per) band is for models where s-wave (p-wave) annihilation 
dominates. 



annihilators and three values of the WIMP spin, ~ 
0, 1/2, 1. We have used the WMAP value of the present 
dark matter abundance, fldmh'^ = 0.112 ±0.009 @. The 
constraint on (Tan is surprisingly robust, as the logarith- 
mic dependence on is offset by the change in the effec- 
tive number of degrees of freedom at different freeze-out 
temperatures. 

Detailed Balancing and WIMP Production at Collid- 
ers — The cross sections of the reaction and its in- 
verse are related by the detailed balancing equation 0: 



= 2 



vii2Sx + If 



(4) 



where the cross sections are averaged over spins but not 
other quantum numbers such as color. For each Standard 
Model particle Xi, we define the "annihilation fraction" 
Ki as 



(5) 



Note that J2 = 1- Using Eqs. Q and ©, we ob- 
tain the following expression for the production of non- 
relativistic x pairs in XiXi collisions 



a{XiXi 

22(Jo-l) 



2x) 
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l/2+,/o 



(6) 



where we have assumed that the inital state particles are 
relativistic {Mx <C M^). This formula is only valid at 
center of mass energies slightly above the 2x threshold, 
V = 2v 



2jl 



4M^/s ^ 1, and receives corrections 



of order v'^. Taking Xi = q ot g (or even W, Z) for a 
hadron collider or Xi = e for an electron-positron ma- 
chine, Eq. ^ provides a prediction of the WIMP pro- 
duction rate. The model-dependence of this prediction is 
contained in a small number of parameters with a clear 
physical meaning: the mass and the spin S^^ of the 
WIMP, the value of Jo, and the annihilation fraction Ki 
for the given initial state. Crucially, the overall scale for 
this prediction, the quantity (Jan, is provided by cosmol- 
ogy, as shown in Fig. 

Tagging and Factorization — Unfortunately, the 2% 
production process whose cross section we predicted is 
not measurable at colliders. Much like Standard Model 
neutrinos, WIMPs cannot be directly observed due to the 
weakness of their interactions with matter. At least one 
detectable particle is required for the event to pass the 
triggers and be recorded on tape. In order to retain the 
model-independence of our analysis, we need to study 
processes in which two WIMPs are produced in associa- 
tion with a photon or a gluon radiated from the known 
initial state. 

In this letter, we concentrate on the case of e+e^ — )■ 
2x + 7- For general kinematics, there is no model- 
independent relation between the rate of this process and 
that of e+e~ 2x predicted by Eq. However, if the 
emitted photon is either soft or collinear with the in- 
coming electron or positron, soft/coUinear factorization 
theorems provide such a relation. Emission of collinear 
photons is described by 



daie^ 



2x + 7) 



dx d cos d 



' J'{x,cos0)a{e+e- 2x), (7) 



where x = 2E.^I ^fs (E^ is the photon energy), Q is the 
angle between the photon direction and the direction of 
the incoming electron beam, T denotes the collinear fac- 
tor: 



Tix, cos ( 



a\ + {\-xY 

TT X 
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(8) 



and (7 is the WIMP pair-production cross section evalu- 
ated at the reduced center of mass energy, s = (\ — x)s. 
Note that upon integration over 0, the above equation 
reproduces the familiar Weizsacker- Williams distribution 
function. The factor T is universal: it does not depend 
on the nature of the (electrically neutral) particles pro- 
duced in association with the photon. Emission of soft 
photons is described by the leading piece of Eq. ||7J) in 
the a; — > limit. 

Combining Eqs. © and O, we find 



da 



dx d cos 9 

aKeiTan 1 -|- (1 - x)^ 



2x + 7) 
1 



167r 



X 



22-^«(2^^-Hl)2 1 



4M. 



2 X l/2+./o 



(1 



(9) 



3 




10 30 30 40 50 60 SO 40 30 BO 100 



E, (GeV) E, (GeV) 



io» 




50 100 150 SOO 50 100 150 SOO 350 300 



E, (GeV) (GeV) 



FIG. 2; Comparison between the photon spectra from the 
process e'^e~ — > 2xi+7 in the exphcit supersymmetric models 
defined in the text (red/dark-gray) and the spectra predicted 
by Eq. @ for a p-annihilator of the corresponding mass and 
Ke (green/light-gray). 
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FIG. 3: The reach of a 500 GeV unpolarized electron-positron 
collider with an integrated luminosity of 500 fb"'^ for the dis- 
covery of p-annihilator WIMPs, as a function of the WIMP 
mass and the e^e~ annihilation fraction Ke. The 3 cr 
(black) and 5 a (green/light-gray) contours are shown. The 
dashed lines include only statistical uncertainty, whereas the 
solid lines include a systematic uncertainty of 0.3% 13- 



This formula, applicable for coUinear photons (6* ^ or 
6' ^ tt), is the main result of this letter. 

Validity of CoUinear Approximation — Eq. ^ pre- 
dicts the rate of events with a single coUinear photon 
and missing energy due to WIMP production in e~^e~ 
collisions. However, very coUinear photons cannot be de- 
tected in an experiment, due to an incomplete electro- 
magnetic calorimeter coverage around the beam pipe as 
well as the lower cut on the Pt,7 = E-y sin 9 that has to be 
imposed to reject backgrounds such as e^e^ e^e^7 
where electron and positron are too forward to be de- 
tected. Do the predictions made using the coUinear fac- 
torization approach have any value in realistic circum- 
stances? To address this question, we have compared the 
event rates obtained by integrating Eq. Q with realistic 
cuts with those obtained in explicit models containing 
WIMPs without making any approximations. For this 
comparison, we have chosen ^/s = 500 GeV, assumed the 
electromagnetic calorimeter acceptance sin6' > 0.1 0, 
and required pr.j > 7.5 GeV corresponding to the mask 
calorimeter acceptance of 1 degree. The results of the 
comparison are shown in Figure |21 The red (dark- 
gray) histograms show the photon spectra from the re- 
action e~^e~ XiXi7 within the minimal supersym- 
metric standard model (MSSM) with the parameters 
suitably chosen to provide the correct neutralino relic 
density [l^. (Explicitly, for = (100,150,200,225) 
GeV, the MSSM parameters take the following values 
at the weak scale: Mi = (115,168,217,242) GeV, = 
(185, 225, 275, 300) GeV, and = (115, 177, 237, 268) 



GeV; for aU four points, M2 = 2Mi, tan/3 = 10, and aU 
the mass parameters not listed above are fixed at 1 TeV.) 
The green (light-gray) lines on the same figure show the 
spectra predicted by Eq. @ for a "generic" p-annihilator 
of the corresponding mass and Ke- We conclude that our 
approach works quite well. The photon spectrum near 
the endpoint is correctly reproduced for all M^. Eq. 
fails for lower values of E^; this effect is especially no- 
ticeable for low M^. This is due not to the failure of 
coUinear approximation, but rather to the fact that the 
relative motion of the produced % particles becomes rel- 
ativistic in this regime, and the higher-order terms in 
the v'^ expansion of Eq. ||2Jl, not captured by CTan, are 
important. Model-independent WIMP searches at e+e" 
colliders, which we discuss below, should take this limi- 
tation into account by concentrating on the photons near 
the endpoint of the spectrum. Note that we did not have 
to impose an additional cut to eliminate central photons: 
coUinear emission naturally dominates the signal. 

Experimental Searches for WIMPs — The main ir- 
reducible background to the search for anomalous 
7-1-missing E events is provided by the Standard Model 
reaction e'^e^ — > lyDj. At the energies well above the Z 
peak, this reaction is dominated by the t-channel W ex- 
change contribution, and has a rather large cross section. 
Nevertheless, the enhancement of the rate predicted by 
Eq. lO may well be observable. In Figure |2| we show 
the reach of a 500 GeV linear collider (LC) with an inte- 
grated luminosity of 500 fb^^ to p-annihilator WIMPs in 
terms of the values of Kg that can be probed at 3 and 5 
a level, as a function of the WIMP mass AI^ . (For com- 
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M^, GeV 75 100 125 150 175 200 225 

cTbg, fb 36 83 202 590 2030 1800 1200 

gBig (p-ann.), fb 1.8 3.9 8.4 21 64 27 4.9 

(Jsig (s-ann.), fb 0.4 0.9 1.9 4.5 13 8.7 3.6 

TABLE L Signal and background cross sections at ^/s — 500 
GeV with no polarization and the cuts specified in the text, for 
a few representative values of M^. The signal cross sections 
are listed for p- and s-annihilators with = 1, and scale 
linearly with Kg. 



parison, a typical value of the Kg parameter in the bulk 
of mSUGRA parameter space is between 0.2 and 0.3.) 
The kinematic acceptance cuts imposed on the photon 
are sin6' > 0.1, pr.^y > 7.5 GeV. Moreover, the accepted 
photons have to satisfy 

The lower cut ensures that the relative motion of the pro- 
duced WIMPs is non-relativistic, < 1/2 in the center 
of mass frame of the WIMP pair. The upper cut cor- 
responds to the endpoint of the photon spectrum for a 
given M-^, and serves to improve the signal/background 
(S/B) ratio. Note that the cuts in Eq. H10|) . and therefore 
the dataset used to test the WIMP hypothesis, depend on 
the assumed WIMP mass. The values of the signal and 
background cross sections with the cuts specified above 
for a few representative values of are given in Tabled 
The table also lists the signal cross section values for the 
case of s-annihilators. It is clear that the sensitivity of 
the proposed search in this case is at best marginal. 

The reach of the LC can be further increased by po- 
larizing the beams: while the background is dominated 
by the W exchange diagrams which only appear for left- 
handed electrons, there is no reason to expect that the 
WIMP couplings have the same asymmetry. For polar- 
ized beams, the signal cannot be fully characterized by 
the spin-averaged annihilation fraction introduced in 
Eq. 0; instead, four independent annihilation fractions 
are needed, corresponding to the four possible e~^e~ he- 
licity configurations. To apply Eq. Q to this case, we 
make a replacement 

«e -> i(l + P-)[(l + P+)^(e^e^) + (l-P+)«(e«e4)] 

+ 1(1 - P_) [(1 + P+He^_e^) + (1 - P+He^^e^)] ,(11) 

where P± are the polarizations of the positron and elec- 
tron beams (P — corresponds to unpolarized beams, 
P_ = 1 to pure right-handed electron beam, and P+ = 1 
to pure left-handed positron beam). Ignoring the Z ex- 
change contribution, the background cross section scales 
as (1 — P-)(l — P+)- For example, let us assume that 



the WIMP couplings to electrons conserve both helic- 
ity and parity, K(e^e^) = n{e^e^) = 2Ke, K(e?e^) — 
K{e^e^) = 0. Then, for electron polarization P_ = 0.8 
and no positron polarization, the S/B ratio is enhanced 
by a factor of 5, whereas for (P_ — 0.8, P+ = 0.6) the 
S/B ratio is enhanced by a factor of 18.5 compared to 
the unpolarized case. 

Discussion — All existing studies discussing the 
prospects for discovery of WIMPs at colliders do so 
within a specific model, fixing the model parameters to 
reflect the observed dark matter abundance. Their inter- 
pretation is hindered by the large number of theoretical 
assumptions made about the details of particle physics 
at the TcV scale. For example, supersymmetric models 
often lead to high rates of events with missing Et at 
the LHC due to production of strongly-coupled gluinos 
and squarks, whose decay chains necessarily involve neu- 
tralinos. It is possible, and indeed likely, that such high 
rates will be observed. However, this is by no means 
guaranteed by the WIMP hypothesis itself. In contrast, 
the rates predicted here are guaranteed (up to the un- 
known annihilation fraction) once the WIMP hypothesis 
is accepted. Therefore, the signatures we have discussed 
provide a unique opportunity to directly test the WIMP 
hypothesis at high-energy colliders. We believe that this 
direction is well worth pursuing. 

Conclusions — We have developed a robust, model- 
independent approach to computing the rates of WIMP 
production in particle collisions. Using the cosmological 
measurement of the present dark matter abundance, we 
have obtained detailed predictions for the rates of pro- 
cesses relevant for WIMP searches at particle colliders. 
Our study highlights the direct connection between cos- 
mology and the physics that may be explored by the next 
generation of collider-based experiments. 
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